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Abstract 


We present an independent catalog (FRIIRGcat) of 45,241 Fanaroff-Riley Type II (FR-II) radio galaxies compiled 
from the Very Large Array Faint Images of the Radio Sky at Twenty-centimeters (FIRST) survey and employed 
the deep learning method. Among them, optical and/or infrared counterparts are identified for 41,425 FR-IIs. This 
catalog spans luminosities 2.63 x 107 < Laa < 6.76 x 10? W Hz ! and redshifts up to z—5.01. The 
spectroscopic classification indicates that there are 1431 low-excitation radio galaxies and 260 high-excitation 
radio galaxies. Among the spectroscopically identified sources, black hole masses are estimated for 4837 FR-IIs, 
which are in 10^? < Mpu € 10?? Mo. Interestingly, this catalog reveals a couple of giant radio galaxies (GRGs), 
which are already in the existing GRG catalog, confirming the efficiency of this FR-II catalog. Furthermore, 284 
new GRGs are unveiled in this new FR-II sample; they have the largest projected sizes ranging from 701 to 
1209 kpc and are located at redshifts 0.31 « z « 2.42. Finally, we explore the distribution of the jet position angle 
and it shows that the faint Images of the FIRST images are significantly affected by the systematic effect (the 
observing beams). The method presented in this work is expected to be applicable to the radio sky surveys that are 
currently being conducted because they have finely refined telescope arrays. On the other hand, we are expecting 
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that further new methods will be dedicated to solving this problem. 
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1. Introduction 


Supermassive black holes (SMBHs, with black hole masses 
ranging 109-10? Mo) are generally identified at the centers of 
massive galaxies with bulges (Fabian 1999; Padovani et al. 
2017). Accretion and ejection are the primary processes in 
which the central SMBHs interplay with their host galaxies 
(Harrison 2017; Zhang et al. 2021), which form the structure of 
accretion disks and jets/outflows simultaneously. Accreting 
SMBHs, a.k.a. Active Galactic Nuclei (AGNs), are highly 
energetic astronomical objects at all electromagnetic wave- 
lengths from radio to gamma ray. In particular, the synchrotron 
emissions from AGN jets dominate the radio band. According 
to the geometrical conditions of jets, radio-loud AGNs can be 
divided into different types, i.e., blazars and radio galaxies. 

Radio galaxies can present compact or extended radio 
morphologies (Miraghaei & Best 2017), with the extended 
galaxies traditionally classified into two main classes, namely 
Fanaroff-Riley (FR) Type I (FR-I) and Type II (FR-ID 
(Fanaroff & Riley 1974). Based on the FR definition, those 
displaying a bright core or an “edge-darkened” characteristic 
along with diffuse lobes are classified as FR-Is, and those 
where the features away from the core are dominated by “edge- 


brightened" characteristics are referred to as FR-IIs. This 
morphological classification scheme was found to be correlated 
with their radio luminosity. Fanaroff and Riley observed that 
almost all the sources in their samples with a radio luminosity 
smaller than 2 x 107 W Hz^! were classified as FR-Is, while 
the brighter sources were predominantly identified as FR-IIs. 
To unveil the optical features of radio galaxies, an alternative 
classification scheme was proposed by Laing et al. (1994) 
based on the intensity of optical emission lines. This scheme 
utilizes the relative intensity of high- and low-excitation lines, 
such as the [O m] emission line, in optical spectra and classifies 
sources into high-excitation radio galaxies (HERGs) and low- 
excitation radio galaxies (LERGs). All reliably classified FR-Is 
were found to be LERGs, while FR-IIs consist of both LERGs 
and HERGs (Buttiglione et al. 2010; Best & Heckman 2012). 
The scheme based on optical spectra can effectively differ- 
entiate between the nuclear and host properties of FR-Is and 
FR-Iis (Capetti et al. 2017a, 2017b). LERGs typically 
correspond to galaxies that experience inefficient accretion of 
hot gas, resulting in the absence of emission lines in their 
spectra. On the other hand, HERGs are commonly associated 
with galaxies that undergo rapid and efficient accretion of cold 
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gas, which is evident from the presence of emission lines in 
their spectra (Best & Heckman 2012; Pracy et al. 2016). 
Therefore, this classification scheme proves beneficial in 
gaining deeper insights into the formation and evolution of 
radio galaxies. 

In the whole FR-II-type radio galaxy sample, those with the 
largest linear sizes (LLSs) are also known as giant radio 
galaxies (GRGs; Machalski et al. 2001; Kuzmicz et al. 2018; 
Tang et al. 2020). They were generally defined as those which 
have LLS close to or larger than 1 Mpc (Willis et al. 1974; 
Planck Collaboration et al. 2016, 2020). GRGs represent a type 
of the oldest radio galaxy and have an essential role in studying 
the AGN feedback process and duty cycle of AGN activities. 
However, the sample of GRGs is relatively small compared to 
the whole population of radio galaxies, which is partially due to 
the weakness of their radio emission. Recently, with the 
operation of high-sensitivity and high-resolution radio surveys, 
such as the LOFAR Two-metre Sky Survey (LoTSS; Shimwell 
et al. 2017), the Rapid ASKAP Continuum Survey (RACS; 
McConnell et al. 2020) and the MeerKAT International GHz 
Tiered Extragalactic Exploration (MIGHTEE; Jarvis et al. 
2016) survey, the GRG samples have increased significantly 
(Dabhade et al. 2020b; Andernach et al. 2021; Delhaize et al. 
2021; Chavan et al. 2023). 

Moreover, straight radio galaxies in FRs are more readily 
defined and computed in terms of their radio position angles 
(RPAs). Analyzing the alignment of RPAs in the orientation of 
radio galaxies across extensive areas of the sky can unveil the 
correlation between the spin axis direction of an SMBH and the 
characteristics of cosmic filaments (Taylor & Jagannathan 2016; 
Contigiani et al. 2017). Osinga et al. (2020) examine the 
uniformity in the position angles of 7555 double-lobed radio 
galaxies (a subset of FR-IIs) using data from the LoTSS 
optical /infrared value-added catalog (Williams et al. 2019) and 
found a significant deviation from the uniformity in samples’ 
two-dimensional (2D) analysis at angular scales of about 4°, 
but no substantial alignment was observed in the samples’ 
three-dimensional (3D) analysis. Simonte et al. (2023) 
investigate the distribution of radio jet position angles in radio 
galaxies within a 1 deg” area of the European Large Area 
Infrared Space Observatory Survey-North 1 (ELAIS N1) field 
observed with the Giant Metrewave Radio Telescope (GMRT). 
The analysis reveals a significant level of spatial coherence at 
position angles exceeding 0°5 and indicates a large-scale 
spatial coherence in terms of angular momentum. 

Given the importance of the FR-II sample, we built a 
machine-learning-based catalog of FR-II radio galaxies from 
the Faint Images of the Radio Sky at Twenty-centimeters 
(FIRST) survey. We performed a cross-match with multi-band 
data, and then physical properties and statistical information 
were derived as well. This paper is organized as follows: 
Section 2 presents the data, search strategy of the candidate of 
FR-IIs and process of finding their optical/infrared 
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counterparts and statistical methods. Section 3 introduces our 
FR-II catalog, provides an analysis and discussion of their radio 
and host properties, explores the GRGs included in our catalog 
and examines the alignment of our FR-IIs. Finally, Section 4 is 
devoted to conclusions. 

Throughout the paper, we adopt a ACDM cosmological 
model with the following parameters: Ho = 67.4 km s~ ‘Mpc’, 
£02,,— 0.315 and Quac = 0.685 (Planck Collaboration et al. 
2020). 


2. Data and Methods 
2.1. The FIRST Survey 


The Very Large Array (VLA) was utilized to conduct two of 
the most ambitious and fruitful radio sky surveys at L-band 
(1.4 GHz): the FIRST (Becker et al. 1995) survey and the 
NRAO VLA Sky Survey (NVSS; Condon et al. 1998). The 
NVSS survey used the VLA D-configuration and achieved an 
angular resolution of 45" and a root mean square (rms) noise 
level of ~0.45mJy, while the FIRST survey adopted 
B-configuration and acquired a resolution of 5" and a typical 
rms noise level of 0.15 mJy. Obviously, the FIRST survey 
surpasses the NVSS survey with nearly 9 times higher 
resolution and 3 times higher sensitivity. Therefore, the FIRST 
survey provides a sharper and deeper view of the radio sky, 
which aids in detailed study of extragalactic radio galaxies. In 
this study, we employed the catalog derived from the data 
release FIRST-14decl7, which contains 946,432 sources 
(Helfand et al. 2015). The FIRST-14decl17 catalog covers a 
total of approximately 10,575 deg? of the radio sky (i.e., ~25% 
of the entire sky), comprising 8444 deg? in the north Galactic 
cap and 2131 deg? in the south Galactic cap. 

The products of the FIRST survey have been extensively 
utilized in the past two decades for studying a diverse range of 
radio galaxies with different morphologies, such as Hybrid 
Morphology Radio Sources (HyMoRS; Gawroński et al. 2006; 
Kumari & Pal 2022), GRGs (KuZmicz et al. 2018; Tang et al. 
2020), winged radio galaxies (Cheung 2007; Cheung et al. 
2009; Yang et al. 2019; Bera et al. 2020), Head-Tail (HT) 
radio galaxies (Missaglia et al. 2019; Pan et al. 2021; Sasmal 
et al. 2022) and double-double radio galaxies (Nandi & 
Saikia 2012). 


2.2. Search Strategy 


The FIRST-14dec17 catalog records the components of radio 
galaxies, and the association and classification of these 
components in radio galaxies need to be done either manually 
or automatically during the post-processing step (Banfield et al. 
2015; Slijepcevic et al. 2022; Lou et al. 2023). We utilize 
HeTu (Lao et al. 2021, 2023), a deep learning-based radio 
source-detector, to automatically locate and classify extended 
radio galaxies within the FIRST images based on their 
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Figure 1. Radio and optical maps of a sample of 12 FR-II radio galaxies. The red color contours represent the radio emission detected in the FIRST image. The 


background color images are from the SDSS r-band. Contour levels are at 30 x [1 


, 1.41, 2, 2.83, 4, 5.66, 8, 11.31, 16, 22.63, 32, 45.25, 64, 90.51, 128, 181.02, 256, 


362.04], where c is the local rms noise. A blue cross is centered on the host galaxy. The background color scheme comes from CUBEHELIX (Green 2011). 


morphology. The latest version of HeTu (HeTu-v2; Lao et al. 
2023) has the capability to classify extended radio galaxies into 
four categories: FR-I, FR-IL HT and Core-Jet (CJ). In this 
morphological classification scheme, we have redefined the 


morphology of FR-II galaxies. We adhere to the original 
definition of FR-II as outlined in the scheme by Fanaroff & 
Riley (1974), but only consider cases where the opening angle 
between the two jets is close to 180°. 
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The FIRST images need to be collected and pre-processed 
prior to source detection by HeTu-v2. To begin with, we 
extracted FITS image cutouts from the FIRST Cutout Server* 
centered on the positions of the sources in the FIRST-14dec17 
catalog. The cutouts measure 3! 96 x 3' 96 (132 x 132 pixels, 
where 1 pixel = 1/8), which ensures that extended sources 
fully contained within fit the cutout and maintain their 
morphological integrity. These cutouts were then converted 
into three-channel (RGB) PNG images using SAOImage DS9 
(Joye & Mandel 2003), utilizing the cool colormap and 
logarithmic min-max scale. Lastly, the HeTu-v2 model was 
applied to both the FITS and PNG formats of FIRST images for 
source detection to build the morphological catalog (called 
FIRST-HeTu), which contains fundamental physical properties 
of detected sources, such as their label, position and flux 
density. A total of 47,582 FR-II galaxies have been detected in 
the FIRST-HeTu catalog. Among them, 45,241 FR-IIs have 
been chosen for inclusion in this study. These galaxies have 
undergone a thorough evaluation and have been shown to have 
a high level of reliability with an average accuracy of 9396 (Lao 
et al. 2023). This paper investigates the further physical 
properties of the selected FR-IIs to build the definitive catalog 
(called FRIIRGCat). 


2.3. Finding Optical/Infrared Counterparts 


We identify host galaxies of FR-II candidates compiled in 
this work, from the Sloan Digital Sky Survey Data Release 16 
(SDSS DR16; Ahumada et al. 2020) and National Aeronautics 
and Space Administration (NASA)/Infrared Processing and 
Analysis Center (IPAC) Extragalactic Database (NED).° FR-II- 
type radio galaxies are those with symmetric jets, and it was 
generally observed that the host galaxies of FR-IIs are located 
at the symmetric center of radio structures (Banfield et al. 
2015). Accordingly, the radio-symmetric centers of FR-IIs 
were cross-matched with the SDSS DR16. Subsequently, we 
utilized the CDS X-MATCH method in TOPCAT (Taylor 
2005) and further identified any potential optical counterparts 
within 15", which is 3 times the angular resolution of the 
FIRST survey. This larger radius was chosen because FR-IIs 
generally exhibit large structures. For any remaining galaxies 
without matches, we selected the source with the smallest 
separation available on NED, using Astroquery (Ginsburg 
et al. 2019) with a search radius of 5”. During these searches, 
multiple counterparts may appear within the designated search 
range. However, following the convention of identifying host 
galaxies for radio sources (Yang et al. 2019), we initially 
eliminated non-galactic objects and then selected the one 
nearest to the central symmetry of FR-II as the final 
counterpart. 


^ https: / [hird.uclInl.org/cgi-bin/firstcutout 
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It is important to note that cross-identifying the host galaxies 
of extended radio sources has always been a challenging issue 
(Alger et al. 2018), particularly for GRGs (Andernach et al. 
2021). Achieving a completely rigorous and accurate determi- 
nation of the host galaxy may even necessitate additional 
optical/infrared observations. Finally, optical/infrared counter- 
parts were identified for 41,425 (38,101 were from the SDSS 
DR16, 3324 were from NED). Nearly 99% of the results from 
NED are infrared counterparts, all of which were derived from 
the Wide-field Infrared Survey Explorer (WISE; Chung et al. 
2011; Rebull et al. 2011; Cutri et al. 2013) data. FR-IIs that 
have SDSS counterparts were investigated for both their radio 
and optical properties. However, the remaining part of the 
study focused solely on examining their radio properties. In 
Figure 1, the FIRST images of a sample of 12 FR-II radio 
galaxies overlaid on the SDSS r-band images are shown, and 
the locations of the optical counterparts of the FR-II galaxies 
are marked with a blue cross (x) for SDSS counterparts. 


2.4. Statistical Methods for Alignment of FR-Ils 


The RPAs for all FR-II galaxies were measured using the 
method described in Osinga et al. (2020). The jet position angle 
is estimated anticlockwisely with regard to the international 
celestial reference frame and the north direction. In this work, 
we define the RPA of jets as the line direction that connects two 
hotspots of FR-IIs, i.e., neglecting the bending effect of jets. To 
assess the deviation from uniformity in the alignment of radio 
galaxies in the sky, parallel transport and dispersion measure 
techniques (Jain et al. 2004; Contigiani et al. 2017; Osinga 
et al. 2020) have been utilized, which take into account the 
effects caused by the geometry of the celestial sphere. First, the 
parallel transport method has been employed to transport the 
"vectors" of the radio galaxy to different locations on the 
celestial sphere. The specific details of this method are 
described below. 

The celestial sphere is constructed using local unit vectors 
(t, Ug, Ug), Where each vector points toward the center of the 
sphere, and northward and eastward on the sphere, respec- 
tively. Figure 2 shows a comparison of the RPAs, PA, and 
PA>, of radio galaxies 1 and 2 with locations, L1 and L2, 
respectively, on the celestial sphere. According to the principle 
illustrated in Figure 2, the comparison of RPAs is converted 
into a comparison of vectors v, and v». 

On a sphere, parallel transport enables us to establish a 
coordinate-invariant inner product between two vectors by 
translating one of them along great circle arcs that connect the 
two (Contigiani et al. 2017). To determine the difference 
between v, and v», the generalized dot product between v; and 
v2 as the dot product between the transported vector vj and v5 
can be defined as 


v On = vi "Wh = cos (PA; = PA> + & = £p. (1) 
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Figure 2. Schematic illustration of parallel transport. This figure illustrates the 
comparison between vectors v, and v». v, is parallel-transported to vj along the 
great circle indicated by the curve from location L; = (7, 01, $4) to location 
L = (2, 92, 2), with u, and u;, being the tangent vectors to the curve at these 
points. The vector v, corresponds to the position angle PA,, while the vector v5 
corresponds to the position angle PA». During parallel transport, the angle 0 
between the tangent vector u, and the vector v, is kept constant. This figure is 
re-plotted from Jain et al. (2004). 


where é; represents the angle between u, and ug, and & 
represents the angle between u,, and Ug,. 

When comparing differences between RPAs, the generalized 
inner product between two RPAs can be redefined as 


(PA), PAs) = cos (2(PA; — PA; + & — &)). (2) 


Considering that RPAs range from 0 to 7, the values of 
(PA;, PA») e[— 1, 1], where +1 expresses perfect alignment 
and — 1 indicates perpendicular orientations of RPAs (Jain et al. 
2004). 

To estimate the significance of a potential alignment in FR-II 
RPAs, the dispersion measure (Jain et al. 2004; Contigiani et al. 
2017; Osinga et al. 2020) will then be employed. The 
dispersion measure of FR-II i as a function of relative RPA 
(PA) is solely determined by the difference between neighbor- 
ing RPAs, which is defined as 


din (PA) = LY pA, PAS, (3) 
R k=i 
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where n represents the number of nearest neighbors considered 
around FR-II i, which includes the FR-II itself. The RPA of 
each neighbor is denoted as PA,. The generalized inner product 
(PA, PA) is defined using Equation (2). 

Then, the non-uniformity of alignment in a sample of N FR- 
IIs is statistically determined by calculating the average of the 
maximum dispersion for a specified number of nearest 
neighbors n, across all N FR-IIs in the sample. This measure 
is expressed as 


1 N 
S, = n di n|max > 4 
DE Z 


where d; nlmax denotes the maximized dispersion around FR-II i. 
Jain et al. (2004) confirmed that in randomly oriented samples 
of radio galaxies, the measure S, follows a normal distribution 
when the conditions N >> n > | are satisfied. This property can 
be utilized to assess the significance of alignment across 
various angular scales. 

Finally, by comparing the statistic S, of the data set with the 
statistical distribution of randomly oriented simulated samples, 
a significance level (SL) is determined to evaluate the rejection 
of the null hypothesis that the galaxy samples exhibit random 
orientation. Therefore, the measurements of S, can be 
transformed into a one-tailed SL when analyzing multiple 
values of n 


SL=1 of 5, (5) 
where 4 represents the cumulative normal distribution 
function, and (S,)uc signifies the expected value for S, in 
the absence of alignment. This expected value is determined 
through Monte Carlo simulations. A lower minimum SL 
indicates stronger evidence for rejecting the null hypothesis, 
suggesting a higher degree of alignment for RPAs. 
Furthermore, if the redshift of the FR-IIs is known, the above 
approach can be extended to identify nearest neighbors in 3D 
space. This enables the detection of potential alignment effects 
and the evaluation of the dependence of S, on a physical scale. 


3. Results and Discussion 
3.1. The New FR-II Radio Galaxy Catalog 


This new FR-II radio galaxy catalog (FRIIRGcat?) contains a 
total of 45,241 candidates. With the available radio and optical 
data, the LLS, RPA, spectral indices, radio luminosities, [O III] 
]uminosities and other host properties of all FR-II candidates 
have been calculated. Based on these results, the statistical 
properties and the uniformity of the alignment of these 
candidates have also been studied. The first 50 rows of FR-II 
candidates in FRIIRGcat are listed in Table 1. The optical/ 
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Table 1 


The First 50 Rows of FR-II Candidates in FRIIRGcat 


Name 


J0000 
J0000 
J0000 
J0000 
J0000 
J0000 
J0000 
J0000 
J0000 
J0000 
J0000 
J0000 
J0000 
J0000 
J0000 
J0000 
J0000 
J0000 
J0000 
J0000 
J0000 
J0000 
J0000 
J000 
J000 
J000 
J000 
J000 
J000 
J000 
J000 
J000 
J000 
J000 
J000 
J000 
J000 
J000 
J000 
J000 
J000 
J000 


+0441 
+0957 
+0958 
—0317 
—0343 
+1002 
—0236 
—0851 
—0808 
+1229 
+0023 
—1121 
+04 
—072 
+143 
—10 

—004 
+12 

—005 
4-072. 
—01 
+02 
+0325 
+0403 
—0939 
+0438 
+0820 
+0545 
+1122 
—0153 
+0324 
+0853 
+0936 
+01 
—020 
+014 
—105 
-11 
—085 
+00 
—04 
407 


SONOS €uro 


[93 
— o R—gERIPO 0) Lk 


RA. 
(12000.0) 


00 00 01.35 
00 00 02.90 
00 00 06.63 
00 00 06.65 
00 00 12.23 
00 00 15.86 
00 00 16.38 
00 00 17.37 
00 00 17.98 
00 00 24.06 
00 00 25.56 
00 00 30.11 
00 00 40.35 
00 00 40.12 
00 00 42.35 
00 00 42.85 
00 00 43.69 
00 00 44.11 
00 00 48.76 
00 00 52.78 
00 00 52.96 
00 00 53.75 
00 00 54.84 


00 0 
00 0 
00 0 
00 0 
00 0 
00 0 
00 0 
00 0 
00 0 
00 0 
00 0 
00 0 
00 0 
00 0 
00 0 
00 0 
00 0 
00 0 
00 0 


02.64 
09.17 
12.10 
14.98 
16.09 
16.58 
18.81 
19.74 
24.81 
26.00 
26.10 
27.28 
30.19 
33.18 
34.15 
34.90 
36.27 
40.86 
43.32 


Decl. 
(J2000.0) 


+04 41 
+09 57 
+09 58 
—03 17 
—03 43 
+10 02 
—02 36 
—08 5 
—08 08 
+12 29 
+00 23 
—112 
+04 1 
—07 2 
+14 35 
—10 13 
—00 4 
+12 16 
—00 50 
4-07 25 
—01 10 
--02 10 
+03 25 
+04 03 
—09 39 
+04 38 
+08 20 
+05 45 
+11 22 
—01 53 
+03 24 
+08 53 
+09 36 
+01 14 
—02 08 
+01 48 
—10 57 
—11 14 
—08 5 
+00 34 
—04 10 
+07 1 


29.6 
00.7 
51.0 
54.8 
01.0 
24.7 
48.4 
23.7 
55.5 
53.2 
28.9 
19.2 
31.7 
32.1 
45.2 
37.1 
30.1 
10.6 
32.7 
23.6 
20.5 
18.6 
36.3 
12.0 
57.1 
55.0 
29.1 
39.3 
21.2 
27.5 
40.6 
25.7 
32.0 
32.9 
41.6 
16.1 
45.4 
12.8 
54.7 
03.7 
30.2 
57.1 


References 


SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
WISEA 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
WISEA 
WISEA 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
WISEA 
SDSS 
SDSS 
SDSS 
SDSS 
WISEA 
SDSS 
SDSS 
SDSS 
SDSS 


RPA 
(deg) 


20.0 
168.1 
26.6 
151.7 
65.2 
141.9 
125.2 
108.5 
39.1 
23.1 
154.5 
147.2 
144.4 
73.3 
161.5 
117.6 
108.4 
26.5 
116.6 
66.8 
16.7 
6.0 
98.7 
100.3 
56.8 
150.2 
30.1 
71.6 
118.1 
120.3 
149.7 
72.6 
23.0 
156.8 
16.7 
116.6 
18.4 
21.0 
114.8 
123.7 
164.7 
47.3 


LAS 
(arcsec) 


21 
35 
24 
27 
26 
32 
37 
23 
37 
69 
42 
30 
31 
132 
23 
47 
23 
28 
40 
14 
19 
34 
24 
40 
62 
29 
129 
17 
61 
25 
25 
30 
78 
41 
19 
16 
125 
25 
26 
26 
21 
32 


LLS 
(kpc) 
143 


225 
242 
236 
197 
270 
497 


57 


161 
338 


96 
62 
205 
157 


846 
109 
328 


250 


82 
208 
171 


Sia 
(mJy) 


22 
306 
306 

40 

60 

71 

63 

23 

41 
109 


S3 
(mJy) 
12 
146 


log(Lraa) 
(W Hz ) 


26.49 


27.68 
27.25 
26.81 
27.40 
27.01 
27.37 


25.33 


26.38 
27.56 


26.36 
25.05 
26.88 
26.63 


27.75 
26.16 
26.30 


27.24 


Dn4000 


1.92 


C, 


2. 


= 


9 
2.31 
Zq 


2.06 
2.23 
29 
3.21 
2.25 
1.62 


2.17 


3.62 


M, 
(mag) 
—22.35 


—25.6 
—21.57 


—26.48 


—22.25 


—2221 


—22.58 


—23.57 
—20.95 


—24.0 


—22.6 


log(Mgu) 
Mo) 


8.41 


8.30 


9.16 


logo m 
(erg s^!) 


39.69 


40.3 


Class 


LERG 


yoy, vpcoc ‘(ddz1) Ic0S€0:vc ‘sorskydonsy pue Áurouonsy ur qoyeosow 
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Table 1 
(Continued) 
Name R.A. Decl. References z RPA LAS LLS Sia S3 apy log(Lraa) Dn4000 €; M, log(Mgu) log(Ljo m) Class 
(J2000.0) (J2000.0) (deg) (arcsec) (kpc) (mJy) (mJy) (W Hz ) (mag) (Ms) (erg s7!) 
J0001—0141 00 01 45.09 —01 41 33.4 SDSS er 157.4 23 er 22 13 0.73 
J0001—1034 00 01 49.93 —10 34 49.0 SDSS er 72.9 24 er 29 16 0.77 er et e er e er e 
J0001--0204 00 01 54.76 +02 04 49.7 SDSS 0.41° 16.4 64 361 309 135 1.09 27.31 e 2.88 —22.28 8.56 e -— 
J0001—0704 00 01 56.01 —07 04 02.0 SDSS 0.44? 157.7 43 253 60 27 1.07 26.67 e er er 
J0001—0944 00 01 59.01 —09 44 58.9 SDSS e 149.0 31 e 19 10 0.84 er 
J0002--1433 00 02 02.51 +14 33 33.9 SDSS 0.49? 50.6 40 250 92 53 0.73 26.91 
J0002.—0002 00 02 02.33 —00 02 44.7 SDSS er 28.8 41 er 22 9 1.16 e 
J0002+1106 00 02 04.99 +11 06 54.4 SDSS 0.40" 90.0 22 122 53 27 0.87 26.49 


Note. Column 1: short name. Column 2: R.A. (J2000, hh:mm:ss). Column 3: Decl. (J2000, dd:mm:ss). Column 4: reference for the position. Column 5: redshift, " represents photometric redshifts, * 
represents spectroscopic redshifts. Column 6: radio position angle. Column 7: largest angular size. Column 8: projected largest linear size. Column 9: integrated flux density at 1.4 GHz. Column 10: 
integrated flux density at 3 GHz. Column 11: spectral index between 1.4 GHz and 3 GHz. Column 12: logarithmic luminosity at 1.4 GHz. Column 13: 4000 A break index. Column 14: concentration 
index. Column 15: r-band absolute magnitude. Column 16: logarithmic black hole mass. Column 17: logarithmic [O m] line luminosity. Column 18: LERG/HERG classification result. 
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Figure 3. Flux densities measured from the FIRST and NVSS are shown for 
5% (2380) of complete FR-II candidates presented in this paper. The red dotted 
line represents the 1:1 ratio between flux densities from the FIRST and NVSS 
surveys. 


infrared counterparts of FR-II sources are listed in the 4th 
column of Table 1, along with the name of the catalog/ 
telescope from which they were identified. The positions of the 
optical/infrared counterparts were used as the positions of FR- 
IIs and are presented in the 2nd and 3rd columns of Table 1. 
When no clear optical counterparts were available, the central 
positions from the FIRST-HeTu catalog were utilized as the 
positions for the FR-II sources, which are indicated as “HeTu.” 

The redshift (z) values of the candidates were estimated from 
the SDSS and NED catalogs and are listed in column 5. The 
redshifts were measured through photometric or spectroscopic 
measurements. For each SDSS counterpart, we first utilized the 
spectroscopic redshift queried from the spectral data of SDSS 
DR16 as the redshift. In cases where the spectroscopic redshift 
was not available, we then adopted the photometric redshift 
queried from the "Photo-z" table in the SDSS DR16 database 
as the redshift value. The redshifts obtained from NED are 
determined by the type code in the “Redshift Flag" column in 
the search results, indicating whether they are spectroscopic or 
photometric redshifts. Redshifts were successfully determined 
for 20,147 out of 45,241 (44.596) FR-IIs, where the number of 
photometric and spectroscopic redshifts were 7261 and 12,886, 
respectively. Most photometric redshifts (12,400) were taken 
from SDSS DRI6. These redshifts were estimated using a 
robust fitting method that considered the nearest neighbors in a 
reference set (Beck et al. 2016). The estimated error for these 
redshifts is at the level of (substantially smaller than) O.1 
(Scranton et al. 2005; Adelman-McCarthy et al. 2007; Beck 
et al. 2016), resulting in large uncertainties in photometric 


Lao et al. 


redshift during 3D alignment analysis (Contigiani et al. 2017). 
Among all candidates, J10534-0604 has the highest redshift 
value, 5.01320 + 0.00167, and J1157+3220 has the lowest 
redshift value, 0.01088 + 0.00001. Four candidates have 
redshift values greater than 4: J10534-0604, J12394-1005, 
J11594-0347 and J15064-4504. The redshifts of these four 
candidates were taken from SDSS DR13 (Albareti et al. 2017). 

Columns 6, 7 and 8 represent the values for the RPA, largest 
angular size (LAS) and LLS of all FR-IIs, respectively. The 
LAS was determined by measuring the distance between the 
endpoints of the two active radio lobes, as measured in the 
FIRST map. In column 9, we chose to measure the radio- 
integrated flux density at 1.4 GHz of FR-IIs using both NVSS 
data and FIRST data. In Figure 3, the distribution of flux 
density for a random selection of FR-II radio galaxies is shown, 
comparing the FIRST and NVSS measurements. This compar- 
ison confirms that a majority (94.6%) of the flux densities 
derived from the NVSS images are notably higher than those 
obtained from the FIRST images. This is because the NVSS 
has more accurate fluxes for extended radio galaxies, allowing 
it to detect low-surface-brightness objects that may be missed 
by the FIRST survey. Additionally, this can be attributed to the 
higher resolution of FIRST and its lack of antennas with short 
spacings. For these FR-IIs, flux density measurements used 
NVSS flux densities. However, a small minority (5.496) of 
NVSS flux densities are smaller than FIRST fluxes. In such 
cases, the flux density measurements used the FIRST fluxes. 
The brightest candidate for FR-II is J12294-0203, with a flux 
(S44) of 55,152 mJy. Conducting higher-resolution observa- 
tions will facilitate the examination of various aspects of this 
candidate, including its dynamics. 

The two-point spectral indexes (oj 4) for all FR-IIs between 
1.4 GHz (NVSS or FIRST) and 3 GHz from the Karl G. Jansky 
VLA Sky Survey (VLASS; Lacy et al. 2020) are listed in 
column 11. The VLASS flux densities listed in column 10 were 
measured from the cutout VLASS images, which were 
obtained by utilizing the cutout core command-line tool’ 
within the Canadian Initiative for Radio Astronomy Data 
Analysis (CIRADA). In FRIIRGcat, spectral index (a; 4) was 
available for 44,763 FR-IIs. The remaining 478 FR-IIs due to 
the VLASS fluxes were not detectable in the CIRADA VLASS 
catalog at that time. The radio luminosity at 1.4 GHz for all FR- 
II candidates is computed in column 12. 

The optical properties of FR-IIs were obtained from the 
SDSS DRI6 spectroscopic catalog and the Max Planck 
Institute for Astrophysics and Johns Hopkins University 
(MPA-JHU, Brinchmann et al. 2004) database. They are listed 
in columns 13 to 17. These properties have been extensively 
described in Kauffmann et al. (2003a, 2003b, 2003c). To 
ensure clarity, we provide a brief introduction to these 
properties as follows: 


7 https://github.com/CIRADA-Tools/cutout. provider. core 
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Spectral index (a3 ,) 


Figure 4. Histogram of the spectral index (a}4) distribution for FR-II 
candidates, showing a mean and median of 0.89 and 0.88, respectively. The 
black dotted curve depicts a Gaussian fit to the histogram. 
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Figure 5. The radio luminosity (Laa) distribution at 1.4 GHz of FR-II 
candidates with redshift (z). The color of each circle in the distribution 
represents the candidate number density with corresponding redshift and radio 
luminosity at 1.4 GHz. The red dashed line refers to z = 4. 


The 4000À break (Dn4000) is an approximation used to 
estimate the age of the stellar population in a galaxy 
(Bruzual 1983). Larger values of Dn4000 indicate galaxies 
with older stellar populations. In this paper, we measure 
the Dn4000 value using the narrow version of the index as 
defined in Balogh et al. (1999). 

Concentration index (C,) is the ratio between the radii that 
encompass 9096 and 5096 of the light in the r-band, which 
can be utilized for the purpose of morphological 
classification of galaxies. Early-type galaxies (ETGs) tend 
to exhibit higher values of C, than late-type galaxies. 
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r-band absolute magnitude (M,) is defined as the apparent 
magnitude in the r-band that the object would have if it 
were viewed at a standard distance of 10 pc (32.6 It-yr). 

Black hole mass (Mpy) is derived from the well-established M— 
c relation, as defined in Tremaine et al. (2002). This 
relation demonstrates a robust correlation between the 
mass of the central black hole (Mgg) and the stellar 
velocity dispersion (0,). 

[O m] luminosity (Lio m) is calculated from the [O m] 5007 
emission line. 


Finally, the HERG/LERG classification results are listed in 
column 18, which used the classification scheme presented in 
Best & Heckman (2012). 


3.2. Radio Properties 


We estimate the power-law (Soc v ^) radio spectral index 
between 1.4 and 3 GHz (oj 4), where S is the integrated flux 
density at a specific frequency v, and a is the spectral index. 
The radio spectral index distribution for all the FR-IIs is shown 
in Figure 4. The FR-II sample in this work has a spectral index 
distribution peaked at 0.89, which is steeper than the typical 
value of 0.7 for radio galaxies (e.g., Kunert-Bajraszewska et al. 
2010). Most of the FR-II radio galaxies (87.0%) show a steep 
radio spectrum with a spectral index greater than 0.5, indicating 
that the measured spectral index is in agreement with the 
theoretical value; 11.4% of FR-IIs exhibit a spectral index 
ranging between 0 and 0.5, which can be attributed to their 
complex structures. The remaining FR-IIs (1.6%) display an 
inverted spectrum and GHz peaked spectrum (GPS) radio 
galaxies, which have small sizes with an average LAS of 23". 

We measured the radio luminosity for all FR-II candidates 
(when both spectral index (o) and redshift (z) are available) 
using the following standard formula in Donoso et al. (2009) 


Lra = AnD?(z)S14(1 + 247}, (6) 


where D, is the luminosity distance to the source in meters (m), 
z is the redshift of the radio galaxy, S;.4 is the measured flux 
density (W m ?Hz !) at 1.4GHz and (1 + z)°~! is the 
standard k-correction used in radio astronomy. 

Figure 5 displays the 1.4 GHz luminosity versus redshift for 
our FR-II candidates with known redshifts. The mean and 
median values of radio luminosity for FR-II candidates with 
available redshifts are found to be log(Lyaq/W Hz ) — 26.61 
and log(Lraa/W Hz !)— 26.57, respectively. This is in agree- 
ment with the division between FR-Is and FR-IIs (Fanaroff & 
Riley 1974). The most radio luminous FR-II candidate is J1239 
+1005, with Lag = 6.76 x 10? W Hz '. Similarly, the least 
radio luminous FR-II candidate is J1109+1206, with 
Lia = 2.63 x 10°? W Hz !. 


Research in Astronomy and Astrophysics, 24:035021 (17pp), 2024 March 


Number 


TM TT Too T 3 1 3T193191311711 


28 


26 —24 -22 


M, 


20 -18 -16 -14 


Lao et al. 


Number 


4 6 8 
log(Mgy) (Mo) 


Figure 6. Distributions of the r-band absolute magnitude M, (left panel) and black hole mass Mgy (right panel). 
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Figure 7. Concentration index C, vs. Dn4000 index for the FR-II candidates (left panel). The absolute r-band magnitude (M,) vs. u — r color for ETG FR-II (right 
panel). The LERG and HERG FR-Ils are represented by green and magenta circles, respectively. The blue dashed line in the right panel represents the threshold that 


separates red and blue ETGs. 


3.3. Host Properties 


The M, values were successfully measured for 7249 FR-II 
candidates. The distribution of M, for these candidates is 
displayed in Figure 6, in the left panel, indicating that the hosts 
cover a range of —14 Z M,Z — 29. Among them, most of the 
candidates (91.8%) have —20 2 M. > — 26. The Mgy were 
successfully measured for 4837 FR-II candidates. The 
distribution of Mgy is shown in Figure 6, in the right panel. 
The majority (92.0%) of them have 9.5 > log(Mgy) Z 7.5 Mo. 

The majority (96.3%) of the FR-II candidates cannot be 
classified spectroscopically due to the absence of emission lines. 
Among the remaining candidates, 1431 FR-IIs are classified as 
LERGs, while 260 FR-IIs are identified as HERGs. Galaxies with 
Dn4000 > 1.7 are defined as red galaxies and galaxies with 
Dn4000 < 1.45 are defined as blue galaxies (Janssen et al. 2012). 
ETGs are defined as galaxies with a more lenient selection 
criterion at C, = 2.6 (Bernardi et al. 2010). Figure 7 shows C, 


10 


versus the Dn4000 index (left panel) for the LERG and HERG 
FR-II candidates. More than 82% of the LERG FR-II hosts are 
located in the region characterized by high values of both C, and 
Dn4000, indicating that they are red ETGs. On the other hand, 
more than 71% of HERG FR-IIs exhibit lower values of Dn4000, 
suggesting that they are blue galaxies. Within this subset, there are 
72 blue ETGs. 

Furthermore, the u — r color of the ETG FR-IIs has been 
considered. Figure 7 (right panel) shows the M, of the hosts 
versus the u — r color. The majority of ETGs are found to lie on 
a tightly clustered red sequence, while a smaller fraction 
exhibits a significant scatter toward bluer u — r colors. We 
applied the method described in Schawinski et al. (2009) to fit a 
Gaussian distribution to the u— r color of the early-type 
population. With the fitted Gaussian, we determined the 3c 
offset to the red sequence. Any ETGs below this offset, 
indicated by the blue dashed line, are defined as blue ETGs. 
Consistent with previous findings based on the Dn4000 index, 
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Figure 8. Radio luminosity at 1.4 GHz vs. [O III] emission line luminosity of the LERG (green circles) and HERG (magenta circles) FR-II candidates. The dashed 
blue line represents the linear fit on HERG FR-II candidates. The black dashed line is the linear fit for LERG FR-IIs with a slope of 0.43. The blue dashed line is the 
linear fit for HERG FR-IIs with a slope of 0.64. 
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Figure 9. Upper panel: histogram of the projected LSS for the FR-II candidates. Lower panel: distribution of the radio luminosity (1444) at 1.4 GHz for FR-II candidate 
LSS. The color of each circle in the distribution represents the candidate number density with corresponding LSS and L,44. The black dashed line represents 700 kpc 
of LLS. 
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Figure 10. Radio and optical maps of a sample of nine GRGs. The contour lines of red, blue and orange colors represent the radio emissions detected in the FIRST, 
VLASS and NVSS images correspondingly. The background images are from the SDSS r-band. Contour levels are at 3c x [1, 1.41, 2, 2.83, 4, 5.66, 8, 11.31, 16, 
22.63, 32, 45.25, 64, 90.51, 128, 181.02, 256, 362.04], where o is the local rms noise. Green crosses are centered on the host galaxies. 


the FR-II HERGs display a bluer color compared to the FR-II 
LERGs. Among the FR-II LERGs, a total of 28 galaxies are 
found to have an u — r color below the dashed line. 

As shown in Figure 8, LERG and HERG FR-IIs exhibit a 
broad range of luminosities in both radio and [Orr] line 
emissions, with a significant concentration around 
log(Lio m/erg s75 = 40. A correlation can be found in LERG 
FR-IIs. The linear fit derived from LERG FR-IIs reveals a 
shallower slope of 0.43, which is less than twice the slope of 
0.99 reported for 3C LERGs in Buttiglione et al. (2010). 
Similarly, a correlation can be observed in HERG FR-IIs. The 
linear fit obtained from HERG FR-IIs indicates a slope of 0.64, 
which is lower than the slope of 1.15 reported for 3C HERGs in 
Buttiglione et al. (2010). 
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3.4. GRGs in the Sample 


The angular size distribution of the FR-II sample in 
FRIIRGcat is from 9” to 214". For those with redshift 
determined, we can estimate their physical sizes. In total, the 
physical sizes (LLSs) of 20,147 FR-IIs are estimated, ranging 
from 4 to 1209 kpc. Figure 9 shows the LLS-L,44 distribution of 
our FR-II candidates. The majority of the candidates share 
modest radio luminosity and LLS. FR-IIs with larger LLS tend 
to have higher radio luminosity. 

GRGs are now generally defined as those with projected 
linear sizes larger than 700 kpc, which yields the classifica- 
tion of 307 GRGs in our FR-II-type radio galaxy sample. We 
note that a few of them have already been listed in the 
existing GRG samples published by Kuźmicz et al. (2018), 
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Figure 11. RPA distribution of the complete sample of selected FR-II galaxies. 
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Figure 12. Illustration of the relationship between the median angular radius 
and the number of nearest neighbors n. The median angular radius is obtained 
by encircling each galaxy with a circular aperture, where the size of the 
aperture is determined by the distance to the nth neighboring galaxy. 


Dabhade et al. (2020a, 2020b), Tang et al. (2020) and Koziet- 
Wierzbowska et al. (2020). We give an overview of the 307 
GRG sample from the literature which results in 5 (J0857 
+3945, J10304-5310, J1236+1034, J1511+0751 and J2356- 
0131), 12 (J1057+5105, J1236--4604, J12564-5528, 71312 
+4809, J13344-5501, J1336+5540, J1341--4916, 71406 
+5547, J1407+5132, J1419+4837, J14264-5129 and 71504 
+5030), 3 (08474-4223, J1129-0121 and J1551+1035) and 
3 (J09414-3126, J1331+2357 and J1421+1016) already 
confirmed GRGs from Kuźmicz et al. (2018), Dabhade 
et al. (2020a, 2020b) and Tang et al. (2020), respectively. The 
identification of GRGs that already exist in historical samples 
confirms the success of our method. The remaining 284 GRG 
sample identified in this work corresponds to the new 
identification. In Figure 10, we illustrate nine GRGs among 
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Figure 13. The distribution of the logarithm of the SL as a function of the 
number of nearest neighbors (n) or the corresponding angular scale is shown 
for the entire sample of selected FR-II candidates. The conversion to angular 
scale is illustrated in Figure 12. 
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Figure 14. An example to highlight the sidelobe residuals enveloping a bright 
source, where the main sidelobes are indicated in white regions. 


the 284 samples, including the largest (J1609+4334) and 
smallest (J15034-2946) ones, which have the largest projected 
sizes of 1209 and 701 kpc, respectively. The new GRG 
sample has redshift span 0.31 « z « 2.42. 


3.5. Alignment in the Orientation of the Samples 


For the sake of maintaining consistency, this particular phase 
of the experiment specifically targeted FR-II candidates in the 
northern sky region, due to observations for half of the southern 
sky region being conducted after transitioning to the new 
Expanded Very Large Array (Helfand et al. 2015). After 
implementing this exclusion criterion, 36,908 FR-II candidates 
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Figure 15. RPA distribution for the sample of 16,828 sources that have a 
redshift measurement. 


were selected, which cover a sky area of ~8444 deg?, resulting 
in a number density of ~4 deg *. 

The distribution of RPAs for all selected FR-II candidates is 
given in Figure 11. Normally, in the absence of systematic effects, 
the RPAs would be uniformly distributed across this large area of 
the sky. In Figure 11, the distribution exhibits four prominent 
peaks around 25°, 65°, 115? and 150°, suggesting the presence of 
systematic effects in these directions. The effect of observing 
beams in the test data comprises the dominant aspect of systematic 
effects. Specifically, the sidelobes in the beam patterns surround- 
ing the bright sources were distributed in these four directions 
(Helfand et al. 2015). Similar effects were observed and discussed 
in Contigiani et al. (2017), which have been proved to be non- 
local, thus having no impact on our analysis. 

To further investigate the conformity of the distribution with 
a uniform distribution of the chosen galaxies, the Kolmogorov— 
Smirnov (K-S) test (e.g., Frommert et al. 2012) was utilized for 
analysis. The K-S test yielded a p-value close to zero percent, 
providing substantial evidence for rejecting the null hypothesis 
that the RPA distribution across the entire selected galaxy 
sample is uniform. 

To estimate the uncertainties and determine the SL in 
rejecting the null hypothesis on different angular scales, we 
compared the results of the selected FR-II candidates with 
those of 1000 simulated position angle samples that were 
randomly distributed. These samples were generated by 
shuffling the position angles among the galaxies while 
maintaining the overall position angle distribution and the 
positions of the galaxies. To represent the statistics in terms of 
angular scale, a circular aperture is plotted with a radius 
extending to the nth neighbor of each galaxy. The number of 
nearest neighbors is then converted to the approximate 
corresponding angular scale by calculating the median angular 
radius of all these apertures. This relationship is illustrated in 
Figure 12. 
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The SL at which the null hypothesis of RPAs being 
uniformly distributed should be rejected is presented as a 
function of the number of nearest neighbors n (or corresp- 
onding angular scale) in Figure 13. There is strong evidence to 
reject the hypothesis of uniformity in RPAs (for alignment) at 
angular scales of approximately 8 deg, with an SL of less than 
10 ^. As evident from the aforementioned RPA analysis, this 
alignment effect is most likely caused by systemic effects 
(Contigiani et al. 2017; Osinga et al. 2020). Here, we present an 
example of the beam pattern in the FIRST image depicted in 
Figure 14, showcasing the sidelobe residuals in a bright source 
due to the sparse uv-coverage of FIRST snapshot observations. 
These sidelobes contribute a significant spurious structure to 
nearby sources and itself (e.g., Helfand et al. 2015). Therefore, 
we attribute the observed alignment in our 2D analysis to the 
sidelobe effect as well. 

Furthermore, the alignment uniformity using 3D galaxy 
positions has also been analyzed. The 3D samples were 
obtained by excluding galaxies lacking redshift information, in 
the northern sky region (36,908 FR-IIs in total). As a result, the 
filer reduces the sample size to 16,828 galaxies. The RPA 
distribution of these samples is illustrated in Figure 15, 
exhibiting similar systematic effects observed in the complete 
sample of selected FR-II candidates. The K-S test also yields a 
p-value close to zero percent, suggesting strong evidence to 
reject the null hypothesis of uniformity of RPAs for this sample 
of 16,828 galaxies. Both 2D and 3D analyses have been used 
for these samples. The statistical method utilized is identical to 
that employed in the complete sample mentioned above, with 
the distinction that the 3D analysis is used to determine the n 
nearest neighbors for each galaxy. 

Figure 16 displays the SL at which the hypothesis of 
uniformity in position angles can be rejected for the 16,828 
galaxies with a redshift. The 2D and 3D analyses indicate 
substantial evidence of alignment at angular scales of 
approximately 075 and 2^0, respectively, with an SL of less 
than 10 ^. The stronger alignment in the 2D analysis, as 
indicated by a lower minimum SL compared to the 3D analysis, 
suggests that the 2D alignment effect is attributable to the 
systematic effects (Osinga et al. 2020). Uncertainties in redshift 
measurements may cause the 3D alignment effect. We 
conducted tests to determine if the results depended on whether 
the redshifts were obtained through photometric or spectro- 
scopic methods. Figure 17 shows the results for the 6132 FR- 
IIs that have a spectroscopic redshift. The figure demonstrates 
that there is no significant evidence of the 3D alignment effect 
in the analysis of these sub-samples. However, in the 2D 
analysis, the alignment effect is still observed. This indicates 
that the presence of the 3D alignment effect in radio galaxies 
can be influenced by less reliable redshift measurements. This 
further substantiates the notion that the 2D alignment signal is 
attributed to systematic effects. These results are consistent 
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Figure 16. The logarithm of the SL at which the rejection of position angle uniformity should occur, as a function of the number of nearest neighbors n, is plotted for 
the 16,826 galaxies with available redshifts. The solid line represents the results of the 2D analysis, while the dotted line signifies the results of the 3D analysis. 
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Figure 17. The logarithm of the SL at which the rejection of position angle uniformity should occur, as a function of the number of nearest neighbors n, is plotted for 
the 6132 FR-IIs with available spectroscopic redshifts. The solid line represents the results of the 2D analysis, while the dotted line signifies the results of the 3D 


analysis. 
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with the findings of the study conducted in Contigiani et al. 
(2017) using the same survey data. 

Several studies (Contigiani et al. 2017; Osinga et al. 2020; 
Panwar et al. 2020; Mandarakas et al. 2021) reported a 2D 
alignment, with angular scales of less than 4 deg. These angular 
scales are similar to those that we have examined in our 3D 
space samples, but half of those we investigated were in the 2D 
space samples. The differences in the 2D alignment angular 
scales were due to the utilization of a different and more 
extensive sample of radio galaxies in this paper. This suggests 
that the 2D alignment of radio galaxies can occur on larger 
scales. However, the reports of a 3D alignment (Contigiani 
et al. 2017; Osinga et al. 2020; Panwar et al. 2020) have not 
yielded statistically significant results. 


4. Conclusions 


In this paper, we present FRIIRGcat, a catalog of FR-II radio 
galaxies consisting of 45,241 candidates from the FIRST 
survey. The code and results for FRIIRGcat are available on 
Github.? A total of 41,425 candidates have been identified with 
optical/infrared counterparts from SDSS DR16 and NED. 
Among them, 20,417 candidates have measured redshifts 
ranging from 0.01 to 5.01. The integrated flux densities at 
1.4 GHz (NVSS or FIRST) and 3 GHz (VLASS) have been 
measured for each FR-II. The spectral index between 1.4 and 
3GHz has a mean of 0.89 and a median of 0.88 for all 
candidates. Most FR-IIs are steep spectrum radio galaxies. The 
candidates with available radio luminosity at 1.4 GHz have 
mean and median values of log(L444) = 26.61 WHz ! and 
logü444) = 26.57 wHz! respectively. 

Most candidates with available r-band absolute magnitude 
have values that fall within the range of —20 € M, S — 26. 
Almost all of the candidates with available black hole mass are 
found within the range of 7.5 S log(Mgy) S 9.5. A total of 
1691 FR-IIs have been successfully classified as LERGs (1431) 
or HERGs (260). More than 82% of LERGs are identified as 
red ETGs, while more than 71% of HERGs are identified as 
blue galaxies. We have discovered a correlation between radio 
and [O III] line luminosity in our HERGs. 

Based on the definition of GRGs, a total of 284 GRGs have 
been identified in our FRIIRGcat. These GRGs cover a redshift 
range of 0.31 « z « 2.42 and were previously not included in 
any existing GRG catalog. These samples will be one of the 
important data sets for the in-depth study of GRGs. 
Additionally, they were automatically identified utilizing the 
results from the deep learning detector HeTu-v2, demonstrat- 
ing the feasibility of employing deep learning techniques for 
the fully automatic identification of GRGs. 

We analyze the alignment in position angles of 36,908 FR-II 
candidates in the northern sky area in FRIIRGcat. The analysis 


5 https: / [github.com/1ao19881213/RGC-Mask-Transfiner/tree / main/tools/ 
inference/FIRST/properties analysis 


16 


Lao et al. 


was conducted on the 2D space (36,908) and 3D space 
(16,828) samples, which were determined based on the absence 
or presence of redshift information in the original samples. We 
discovered evidence of alignment in both 2D and 3D space 
samples. The 2D alignment is most likely due to systematic 
effects (e.g., the strong sidelobes in beam patterns), whereas the 
3D alignment is likely caused by the less reliable redshift (i.e., 
photometric redshifts) because there is no alignment signal that 
was observed in the 3D sample with spectroscopic redshifts. 
The angular scale of the 2D alignment in this work is twice that 
of the previous works. In the case of reliable redshift 
measurements, this study confirms that the alignment of radio 
galaxies identified from the FIRST survey is primarily 
attributed to systematic effects. This work provides important 
insights for further investigating the cosmic alignment of AGN 
jets, which are indicators of large-scale structures and important 
for cosmological models and parameters. These effects will be 
important for cosmological analyses involving radio data, such 
as weak lensing studies conducted with advanced high- 
resolution radio interferometers like the Square Kilometre 
Array (SKA; Weltman et al. 2020). 
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